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I. INTRODUCTION 

The investigation of energy spectrum of light muonic atoms (muonic hydrogen, muonic 
deuterium, ions of muonic helium) reached a new level at present. This is due to new experi¬ 
mental results obtained by the CREMA collaboration in [l|-l3| . On the one side these results 
open a possibility to obtain new values of a number of fundamental physical constants such 
as nuclear charge radii. While experimental data on atomic transitions have become very 
precise, our knowledge of the charge radii, which are part of theoretical predictions, is not 
as accurate as we would like. On the other hand, they call to look again at the formulation 
of the theory of bound states in quantum electrodynamics and possibly revise some of its 
previous aspects. The second position was proved important after a series of experiments in 
[I|, which revealed essential disagreement between two values of the proton charge radius 
obtained in experiments with electronic and muonic atoms 0,1,1. An analysis of the situa¬ 
tion and determining the causes of discrepancies are investigated in several directions, which 
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are widely discussed in j3, |5-13|. It is possible that the publication of new experimental data 
on the structure of the energy levels of muonic deuterium which is planned in near future, 
will help clarify the problem. A comparison of the theory and experiment for the transition 
frequencies 2‘^^'^^Sji) in muonic deuterium demands careful consideration of 

different contributions to the energy P-levels. The calculations of hne and hyperhne struc¬ 
ture of the energy spectrum of light muonic atoms were made in a series of papers [li-H 


The results of these studies are a reliable benchmark for a comparison with experimental 
data and provide a starting point for further research. Whereas the calculation of separate 
contributions to the hyperhne structure of S-states of the muonic deuterium, even with a 
very specihc kind, was the subject of intense study, the hyperhne structure of P-states much 
less investigated. Therefore, in this study we aim to partly hll this gap. In this work we make 
new analysis of diherent corrections to hyperhne splittings of P-states which allow to obtain 
more accurate results important for a comparison with experimental data. Another aim of 
our study is to develop a method of projection operators in the investigation of the energy 
structure of P-states. The method of projection operators on the bound states with dehnite 
spins was used previously in [l3, [isj for the construction of particle interaction operator for 
hyperhne structure of S-states. 


II. GENERAL FORMALISM 


Let us begin our consideration with basic contributions to hyperhne structure of P-states 
of order Our approach to the calculation of hyperhne splittings is based on quasipotential 
method in quantum electrodynamics in which the two-particle bound state is described 
by the Schrodinger equation [liQ. In this work we develop another approach to the 
calculation of hyperhne structure of muonic deuterium based on tensor representation of P- 
wave projection operators describing muonic deuterium states. First we show on an example 
of calculating the leading order contributions how a tensor formalism helps investigate the 
hyperhne structure of the spectrum. It is useful to work in momentum representation where 
we can write the wave function of muonic deuterium 2P-state in the tensor form: 


V’2p(p) = {s-np) i?2i(p), 


( 1 ) 


where es is the polarization vector of orbital motion, Up = (0, p/p), i? 2 i(p) is the radial wave 
function in momentum space. Then the energy shifts are presented in integral form: 

= / (e- . n,) j (e ■ „,) q). (2) 

In the leading order the hyperhne potential is constructed by means of one-photon 

interaction amplitude Ti^ . Writing the amplitude Ti^ we refer to it a part of the bound 
state wave function related to orbital motion: 


^17 (P) q) = dvrZa (e* • Uq) 


“(^ 1 ) + (1 + “(Pi) 


* ( \ f (P2 + ^ 2)1 


2mi 

(P2 + ^ 2 ) 1 / kpk 


2mi 


ky 


{e • Up) Dpuik^ ^ 

(3) 


Fi{k^) - ^ ^ ^ -"A 17 ^n2 


F2{k ) T (PpAPcrp 9pli9u\) F^{k ) |■£'c^,c^(P2)^ 


2m2 2m\ ^ 2m 
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where pi ^2 = (m 7 +m 2 ) ~^^^ four-momenta of initial muon and deuteron, gi 2 = 
q are four-momenta of final muon and deuteron. They are expressed in terms of total two- 


particle momenta P, Q and relative momenta p, q. is the photon propagator which is 

taken to be in the Coulomb gauge. Explicit expression of the deuteron wave function ed,a{p) 


has the form: 


^dAp^) 


1^2,(7 + goam 2 {eaA^) ' P 2 ) 

/ N 

e2[p)+m2 m 2 


(4) 


It should be noted that the amplitude ([3]) has been studied in detail in |22| excepting 
quadrupole correction. In the center-of-mass rest frame P = Q = Mv, v = (1,0). The form 
factors El 2 , 3 (/c^) are related to the charge, magnetic and quadrupole deuteron form factors 
as (p = A/4ml) in,Hi]: 


Fc = Fi + ‘^v[Fi +{^ + v)F 2 -Ps], Pm = F 3 , Pq = Pi + {I + p)P 2 - Ps- (5) 

We consider ((3]) as a starting point for a composition of orbital L momentum, the deuteron 
spin S2 (note that the spin of the nucleus is usually denoted by I) and muon spin Si. In 
the hrst scheme of momentum composition we add hrstly momenta L and Si obtaining two 
muon states with angular momenta j = 1/2 and j = 3/2. In the Rarita-Schwinger formalism 
the wave function of the state with half-integer spin 3/2 is described by 

^/.(p, A)mc^(p,w), (6) 

where (la;;lA | |(j) are the Clebsch-Gordon coefficients. Another sequence of angular mo¬ 
mentum addition is that in the beginning we add the orbital and intrinsic angular momentum 
of the deuteron and then the muon spin. When we combine the L = 1 and S 2 = 1 we get 
three states with the deuteron momenta 2,1, 0. The deuteron wave function has in this case 
the form: 

7) = 57 ^^2 I 27) eAp, Ai)£^(p, Ai). ( 7 ) 

Ai,A2 

After combining with the muon spin on the second stage there arise three states with 
F = 5/2, 3/2,1/2 which are described by the tensor-spinor held 4/^^ satisfying to the Dirac 
equation 

(h - 1)T^^ = 0, = 0. (8) 

The held can be easily decomposed into diherent parts with dehnite atomic angular 
momentum F\ 

^’'■ = 1 : 4'^.; (9) 

= , ( 10 ) 

= \ '■ ^ , (11) 
= \ ■ = ^ 7±^], 



( 12 ) 
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F^ = 




flU 


Vs 


{g^,u - Vf,v^), 


(13) 


where is the usual 5/2 generalized, symmetric Rarita-Schwinger tensor-spinor [251. l26 


The negative parity is obvious from physical reasons. Different states with total momentum 
F = 1/2 and F = 3/2 are decomposed into symmetric and antisymmetric parts satisfying 


to dS]). The tensor-spinor wave functions were used previously in j27|] for the bound states 
of quarks. For further calculations, we note that each held with F = 3/2,1/2 is a 
superposition of states with muon angular momentum j = 1/2 and j = 3/2. Introduced in 
tensor-spinor helds can be considered as both a projection operators on the states 
with a dehnite value of the total angular momentum. These projectors are very convenient 
for the calculation of the matrix elements of the interaction potential corresponding to certain 
quantum numbers. They allow us to avoid direct cumbersome multiplication of different 
factors in the amplitudes of the interaction of particles and use the computer methods for 
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calculating amplitudes and the energy shifts 

To demonstrate this property of we continue our calculations of the amplitude ([3]) 
corresponding to transitions between states with dehnite values of F. Introducing projectors 
in dSD and averaging the amplitude over the projection of the total angular momentum 
Ai we obtain the following basic relation: 


Ti^{p,V) = 


- n„n„lr< 

2F + 1 ^ ^ 


E 

M=-F 


UJ(7l dpi 


[(mi(h + 1) - 7q]^ [mi{v + 1) --fp] 


2mi 


2mi 


X 


9po 


iP2 + 92)1 


2m2 


FiiV) 


{P2 + 92 )u kpK 
2ml 


2m2 


F2{V) + {Ppxg, 


ap 


gppgax)7^F3{V 

2mo 


(14) 




gppi 


21712 


{m2Vp^ - qpP{2m2Vp - qp) 


gaai 


27712 


{m2Va^ -PaP{2m2Va - p^) 


where the lepton vertex function T^ = + (1 + cip is the muon anomalous 

magnetic moment. The Lorentz factors of the Dirac bispinors and transformed Lorentz 
factors of deuteron wave functions are written explicitly. Inserting in ffTTj) from (p|)- 
m, averaging and summing over initial and hnal state polarizations Ai and calculating 
the trace by means of the package Form 2^ we hnd three matrix elements corresponding 
to F = F = I F = The polarization sums for the helds with half-integer spin looks as 
follows 
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n^.(F = 3/2)= 

M=-F 


M 

V 


(h + 1) 


1 2 1 , 
gpv - IpiAv - - Vu-ip) 


, (15) 


M=-F 


M 

pa 


(h + 1) 


_( pi + F^ F^ 1 - -F^ F^ - 

2 y pp Fa ' pa Fp j 2 pF pa 


( 16 ) 


_+ P^P^P^ + P^P^P^ ) P^ = a 

10 p ^ va ' ^ V ^ p ^ pa ' ^ a^ vp ' ^ u ^ a^ pp) i ^ pu iipp 


VpVy, 


F^ = F^ 0/ 

p pF /M* 


Let us construct by this method basic hyperhne splittings of order a^. We project the 
amplitude ([3]) sequentially on states with j = 1/2, F = 1/2 and j = 1/2, F = 3/2. 
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Corresponding averaged amplitudes are the following: 


rF=i/2 irZa 5 


7^17 (p> q)j=i /2 = + l)( 7 pi - v )(75 + Vs) 


9 


[{mi{v + 1) - 7q] 


2mi 


X 


(17) 


xr, 


[mi{v + 1) - 7P] 

2mi 


( 7 ^ + - 'y<7i) \D^^{k)x 


k> 


^ 3pf^9ax)^^F^{k‘^)\x 


X 


5'ppi 


1 


2?77n 


{rn2Vp^ - qp^){2m2Vp - qp) 


1 


dual 


2m: 


im2Vai -pai){2m2Va -pa) 


_ ^ 3/2 TT Z ex r 

Ti^ip, q)j=i/2 = + l)naipi(i^ = 3/2)(75 - n5)75X 

^ [(m.(f. + 1) - 7qi r Jm.(f.+ 1) - 


(18) 


27771 


2m^ 


9 pi 


k> 




X 


9ppi 


27772 

1 

27?7n 


{m2Vp^ - qp^){2m2Vp - qp) 


1 


Qucri 


2m: 


im2Vai -pai){2m2V^ -pa) 


In the quasipotential method each of the amplitudes (ITTll - ffTSl) determines the interaction 
operator of particles corresponding to states with selected quantum numbers. In this case, 
we get not only the contributions of the hyperhne interaction, but also the Coulomb potential 
and a potential of hne structure. But the difference flTTl) and flTSl) allows to hnd the hyperhne 
splitting of state j = 1/2 which is written as an expression of the output from the Form 
program: 


zhfs 


Za { 


rn(P,q),:i/ 2 (P = 3/2;l/2) = — 


+ (^d + 1) 


2 (pq )2 (p^ + g2)(pq) 




7771 


\rn 2 Kd 


'_m (pq)^ ^ 


+ 


(19) 






+ 2il + Kd)il + J) 


M 


(pq) 

pgk^ 


+ 


+ g^)(pq) 

pgk^ 


where Kd = 0.714025pAr is the deuteron anomalous magnetic moment [29|, connected with 
the deuteron magnetic moment pd by the relation Hd = {pdm 2 /mp — 1). In flT^ we take 
electromagnetic form factors at = 0 and omit the quadrupole contribution which is 
studied in detail in next section. Normalization factor S/dvr coming from wave function of 
orbital motion is taken into account. Two other hyperhne splittings of 2 P 3/2 state looks as 
follows: 


-hfs , , , X Za 777i 

'iUpp = 3/2; 1/2) = ^ ' 

(pq)^ l(pi + 9^)(pq) 


7’n(p,q),:3,2(7 = 3/2il/2) = 

+(x;d+l) 


pgk^ 


pgk^ 


7772 Kd 
2 


Ipg ^ 11 (pq)2 3 (p2 + g2)(pq) 


2 k^ 10 pgk^ 


10 


pqk'^ 


+ ( 20 ) 


+ 5(l+xj(l-^) 


M 

■k2 


(pq)^ ^ (p^ + ?^)(pq) 


pgk2 


pgk2 
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+ (^(i+l) 



-(1 + ^d) — 

10 7772 

(1 -|- Kd) ( 

pq 

[2 

5 v 


7771 

^ 7772Krf 

'5pq l(pq) 
6 k^ 2 pgk^ 

5(pq)2 5(p2 + g2)('pq)- 

3 pgk^ 6 pgk^ 

+ |(1+kJ(i-y) 


I(p 2 + g 2 )(pq) 


+ ( 21 ) 


(pq) 

pq 


El + 

k2 pgk2 




pgk2 

(pq)2 (p2 + g2)(pq) 


pgk2 


where the terms proportional to {pc[/pq) vanish as a result of the angular integration. They 
are important for the correct calculation of the vacuum polarization effects. There are three 
types of integrals with the radial wave functions which are calculated analytically: 


Ji = 


^21 (g) 


dq 

(27r)3/2 


^21 (P) 


dp pq 

(27r)3/2 (p _ q)2 


PQ \ ^ 

(p-q)v 16’ 


( 22 ) 


(pq)^ \ 

pg(p - q)v 


^ / (pq)(p^ + g^) \ 

48’ ^ \ pg(p-q)2 / 


5 

2A' 


Note, that the terms on the right side of the equations (lT^ - (|2lil proportional to (I + k^) dis¬ 
appear after momentum integration and we obtain the following leading order contributions 
of diagonal matrix elements to hyperhne splittings of 2 Pi /2 and 2 P 3/2 states: 


= 3/2; 1/2) 

= 3/2; 1/2) 
A£:‘ 42 (f’ = 5/2;3/2) 


0)^(1 + Kd)fi^ ^ ^ miKd ^ 

12 mim 2 2777,2(1 -|- Kd ) 2 


a;'^(l Kd)p^ 

'2 

5 + 

miKd 


24777i7772 

27772(1 Kd) 

10_ 

5q;^(1 Kd)p^ 

'2 

.5 + 

miKd 


7277717772 

27772(1 -F Kd) 

10 


2070.5040 peV, 

■ 420.9426 peV, 
= 701.5712 peV. 


(23) 

(24) 

(25) 


Two amplitudes ffT7)l - ffT8|) are constructed combining hrstly the orbital momentum and 
muon spin. Then the spin of the nucleus is added. We can act slightly different express¬ 
ing the states with j = 1/2 and j = 3/2 directly in terms of introduced symmetrical and 
antisymmetrical states. This possibility is illustrated hereinafter. In this method we can 
evaluate also off-diagonal matrix elements. Their calculation is demonstrated in next section 
for quadrupole correction. To facilitate a comparison of the method of calculation and ob¬ 
tained contributions to the previous approaches we make Appendix A, which demonstrates 
the calculation of corrections of order ck^ in the coordinate representation. All basic con¬ 
tributions to hyperhne structure and numerous higher order corrections are presented in 
Tabled! 
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TABLE I: Diagonal matrix elements of hyperfine structure of 2P-states in muonic deuterium 


The contribution 

2^Fi/ 2, peP 

2^Fi/ 2, peP 

2^F3/2) peP 

2^F3/2) peP 

2®F3/2) peP 

leading order 

correction 

-1380.3360 

690.1680 

8162.2889 

8583.2315 

9284.8027 

quadrupole correction 

of order 

0 

0 

433.9033 

-347.1227 

86.7807 

vacuum polarization 

correction of order 

-1.0706 

0.5353 

-0.2802 

-0.1121 

0.1681 

quadrupole and vacuum 
polarization correction 

of order 

0 

0 

0.3564 

-0.2851 

0.0713 

relativistic 

correction of order a® 

-0.1677 

0.0838 

-0.0125 

-0.0050 

0.0075 

vacuum polarization 
correction of order a® 

-0.0011 

0.0005 

-0.0014 

-0.0006 

0.0008 

structure correction 

of order a® 

-0.0011 

0.0021 

-0.0006 

0.0010 

-0.0016 

Summary contribution 

-1381.5765 

690.7897 

8596.2539 

8235.7070 

9371.8295 


III. QUADRUPOLE INTERACTION CORRECTIONS 


Quadrupole interaction originates from not completely spherical shape of the deuteron. 
If the potential of the muon has also a non-spherical component at the position of the 
deuteron (only with muon angular momentum j > 1/2) then there exists a quadrupole 


energy shift |3 Ch 32| . Ordinary calculation of this contribution to hyperhne structure in 


muonic deuterium is based on the representation of quadrupole interaction in coordinate 
space as a scalar product of two irreducible tensor operators of rank 2. After that the 
matrix elements of tensor operators are expressed in terms of reduced matrix elements using 
the Wigner-Eckart theorem. 

In this work we develop another approach to the calculation of quadrupole interaction 
based on tensor representation of P-wave projection operators describing muonic deuterium 
states. In the case of F = ^ and F = | we should take the sum of two contributions regarding 
to symmetric and antisymmetric wave projection function fll0p - fll3p . For completeness, we 
present two averaged amplitudes corresponding to = |) and = |): 


TTa 


Tiy{p, q)^‘ = (h + l)n^5(7p - Vp) 


X 


x^9pc 


4 ■■'q ~p 

[mi(l -v) + 7 p] 

2mi 

{P2 + Q2)u t:, n 2 


[(mi(l -v) + 7 q] 

2mi 


r^x 


(26) 


(Tcj - ^'a;) + (p -t 5, w -)■ a) - (ca ^ cr) - (p -)■ d) IX 


























pq 


ipqj 


E A- 1 
q p 


(p - q)^ 


(p - q)^ 


+ 


1 (pq) 


(pq)^ _ 

p - q)2 3 pq 


Ti^{p, q)^ = ^nlnpDf,^{k)Tr\ {v + l)n^5(7p - Vp) 


[(mi(l -v) + 7 q] 


2mi 


xigpc 


[mi(l --0) + 7p] , N / e- \ / N / 

X--( 7 ij — Vuj) + {p ^ 6, u ^ a) + {u a) + {p 

{p2 + q2)u {P2 + q2)ukpK 


r^x 


X 


(27) 


2m2 


Fi{k^ 


2-1712 2mi 


m") + (SmS.;. - S»9.a)t^^F3(*:") } = 


2m2 


7i(yC^d 

15 


pq 


(p - q)^ 


(pq) + § 

(p - q)^ 


(pq)' 


1 (pq) 


(P - q)^ 3 pq 


where we keep only the contribution of the quadrupole form factor Fq(0) = Qd. The 
index replacements designated in brackets of fl26l) and fl271) refer to the written part of the 
amplitude. Remaining integration with the radial wave functions is carried out analytically: 


J = 


dp 


(27r)3/2 


Rl2{p) 


dq 


(27r)3/2 


Ri2{q) 


pq 


(p - q)^ 


(pq) U ; 

(p - q)^ 


+ 


(pq)^ 

(p - q)' 


p^{Zay 


Ibvr 


(28) 

The sum of fl26|l and flTT)) multiplied by the factor fl28l) gives the contribution to hyperhne 
splitting pj-gggnt final results for other transitions: 


ABy* = 


di^pZ cr) 
48 


5^1 ~5p3 + —Sp5 

2 5 2 5 ^2 


(29) 


The quadrupole moment of the deuteron is taken to be Qd = 0.285783(30) /m^ [3^. The 
result fl2^ coincides exactly with previous calculations made by different approaches fl^ . 
As it follows from numerical values of fl29|) (see Table [T|) the quadrupole interaction changes 
the position of levels 2^P3/2 and 2‘^P^i2. Let us investigate in addition how the total angular 
momentum of the muon is changed in such transitions. For this purpose, build again the 
amplitude of single-photon exchange combining consistently muon spin with the orbital 
angular momentum and the deuteron spin. To be specihc, we consider two diagonal matrix 
elements which are determined by averaged amplitudes with j = l/2, P = l/2 and j = 1/2, 
F = 3/2: 

Ti7(p,q)ii" =—nlnpDp^{k)Tr^{y + l){pip-Vp){pis + Vs)x (30) 


X 


[(mi(l + h)- 7 q] [mi(l + h)- 7 p] 


2mi 


2mi 


x\gpc 


{P2 + q2)v 
21712 


{P2 + q2)u kpk. 


{'lL2+Vuj)ha-Va)\x 


k\ 




rF=S 


ttq; 


rii(p,q)ii^ ^ —nji'tDi^{k)Tr\ (v + TjU^piys - vs)25'-< 


(31) 
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X 


[(mi(l + {})-7q]T. [mi(l+{)) - 7p] 


X 


9po 


{P2 + 92)1 


■F^{k^) - 


2mi 
{P2 


2mi 


75 ( 7 ^ - [x 


(I 2 ) 


h h 

V 


k\ 


F2{,k) + {gpxga^, - gppgax)^. — F^{k ) ^ = 0 , 


2m2 2m2 2771^ p p pp j 2m2 

where lower indexes of the amplitude designate the muon total angular momentum. With 
one side, the obtained expressions fl5UD - fl?Il) explicitly show that quadrupole interaction 
does not contribute to diagonal matrix elements with j = 1/2. On the other side they 
demonstrate our choice of the tensor projectors on the state with j = 1/2: 


' (j = 1/2) = ^75(7 m - w)75(7!^ - (32) 

where the spinor 'h describes the state with total atomic momentum F = |. Using the Dirac 
algebra transformations we can expand fl3^ on the basis = ^) and = 2); 


^(j = 1/2) = = 1/2) + = 1/2). (33) 

The same expansion can be performed for the state with j = 3/2 and two states with J = |, 
F = I and j = F = |. They looks as follows: 


>(/ = 3/2) 



= 1 / 2 ) 



U(F = 1/2), 


(34) 


U. Hj = 1/2) = = 3/2) - yfU/i" = 3/2), (35) 

<„■"(/ = 3/2) = = 3/2) + ]IIK<.P = 3/2). (36) 

Using fl3^ - (l36|) we can investigate off-diagonal matrix elements corresponding to different 
values of muon angular momentum j. In fact, contributions with symmetric and antisym¬ 
metric tensor-spinor fields |) and = i, |) are evaluated above in matrix 

elements fl26|) - fl27)) . Thus it is necessary to use only the correct coefficients of expansions 
fl33|) - (l36|) . As a result we obtain: 

AE»‘U = 1/2;j' = 3/2) = (y2S^, - . (37) 

Numerically, all quadrupole corrections are large and presented in Table [I] and Table [TTl 
Drawing attention to the signihcant value of the quadrupole corrections, we proceed to the 
consideration of other important effects within the formulated framework. 


IV. VACUUM POLARIZATION AND STRUCTURE CORRECTIONS 

The above basic formulas for the amplitudes of the muon-deuteron interaction allow to 
calculate the various corrections. Next in importance are the corrections to the vacuum po¬ 
larization (VP) of order a^. In the formulated framework these effects can be easily studied. 
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In the first order perturbation theory one-loop vacuum polarization contribution to HFS 
is determined by the amplitude in Fig. [1] For its calculation in momentum representation 
which we use, the following replacement in the photon propagator should be done in ([T9]): 


1 « r Pirn 

k? Stt k? -|- ’ 


p(?) = \/F^(2y+i)/y. 


(38) 


As a result we hnd that the vacuum polarization contribution to hyperhne splittings can 
be expressed in terms of three momentum integrals which are a generalization of the three 
integrals discussed earlier in fl22|) : 


h 


/ 


R2i{q) 


dq 

(27r)3/2 


^21 (P) 


dp pq 

(27r)^/2 (^p _ q^2 _j_ 4 j7^2^2 


(39) 


pq 


h = 


(pq)" 


pg(p — q)2 -|- 4m2,^2 


a(3a -|- 8) -|- 6 
(p — q )2 -I- 4 m 2,^2 y/ 2(a-|-2)"^ 

a(3a -|- 8) -|- 10 ^ / 


a = 




6(a + 2)4 


h = 


pa 

(pq)(p2 + g2) 


\pq{p — q)2 -|- 4m2,^2 


2(4a + 5) 
3(a + 2)4' 


Third integration over the spectral parameter ^ also can be carried out analytically, but they 
are quite cumbersome. So, we present here necessary VP correction to hyperhne splitting 
of 2 Pi /2 state only in integral form: 


AE»’(2P,/2) 


p?a{ZaY 

37rmim2 



'rriiKd (3a -|- 2) 

- 2m2 3(a -|- 2)^ 


h(l + /t(i) 



~(1 + ^d) 


2(a + 2)4j 


1.0718 peV. 


2(3a2 + 4a + 2) 
3(a + 2)4 

(40) 


The same calculation can be performed for the 2 P 3/2 state. The corresponding results are 
the following: 


AE;P(2P3,2)(F = 3/2;l/2) 


p'a{Za)^ P (3a + 2) 

67rmim2 ^ - 2 m 2 6(a -|- 2)^ 


(41) 


, ^ P (15a2 + 8a + 4) 

+ (1 + K.) (I-tJ 3(a + 2)* 


(1 + ^d) 


2a^ 1 

(a + 2)V 


0.0595 peV, 


A-B;/'(2P3/2)(F = 5/2;3/2) 


p^a{ZaY f°° VTTiiKd 5{3a + 2) 
67117111712 Ji ^ - 21712 18(a -I- 2)3 


(42) 


+ n + K 1 ("l - 2(15a^ + 8 a + 4) 

4 j 9(a + 2)4 


(1 -|- Kd) 


lOa^ 

3(a + 2)P 


0.0992 peV. 


Another important VP effect is related to quadrupole interaction discussed in previous 
section. Using for its calculation basic expression flT^ . fl20|) . (1^ and (l39|) we obtain for 
diagonal and off-diagonal matrix elements: 




p^a{ZaYQd f°°{5a^ + Sa + 4:) 


367r 


(a+ 2 )^ 




4 

5^ 


Srpl 17 5 

J- 2 " o 


F- 
^ 2 


( 43 ) 
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AEq,,p{j = 3/2; j' = 1/2) 




1 

— ( 
5 


+ T^F- 
2 


X 0.2441 /ieV, 


fi^a{ZaYQd P (5a2 + 24a + 24) 

-A 3(a + 2)- 




1 

7 ! 



(44) 


FIG. 1: Vacuum polarization effects in one-photon interaction. The wavy line represents hyperfine 
part of the interaction. 


A comparison of our results fj4ni) - fj4411 with earlier estimates in shows that there is a 
significant difference of the order of tenths of fieV. For this reason we decided to perform 
additional validation of our results using a different method of the calculation. As was shown 
in jlll the vacuum polarization effects presented in Fig. [T] in first order perturbation theory, 
can be calculated in coordinate representation. The amplitude shown in Fig. [T](a) gives the 
following hyperhne interaction potential in coordinate space: 




Za{l + Kd) « 




—2meir 


1 + 


miKd 


2mim2r^ Stt m2{l + Hd), 

X (L ■ S2)(l + 2me^r) - {I + a^) ■ S 2 ) - (si ■ n)(s2 ■ n)]-F (45) 


+ {1 2me^r)[{si ■ S 2 ) - 3(si ■ n)(s2 ■ «-)] 


Averaging fl45|) over the Coulomb wave functions, we obtain an analytical expression for the 
vacuum polarization correction of order in one-photon interaction: 




24mim2r'^ Ptt Jq 


1 + 


miKd 


m2(l + Kd) 


X 


—, 2mnj\eP . , . f — , 2me£ 

xT.d + - (1 + a,) + (1 + j 


where we introduce the designations for operators Tj in 


T, = {L-S2), T2 = 


{si- S 2 ) -3{si-n){s2-n) 


n = 


[si- S 2 ) - {si-n){s 2 -n) 


(46) 


(47) 
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The coordinate integration in (H6|) is carried out analytically and numerically over the 
spectral parameter Numerical results for separate states include both diagonal and off- 
diagonal matrix elements: 

= 1/2) = -0.7145 (48) 

= 3/2) = 0.3573 fieK 
= 1/2) = -0.0992 ^,eV, 

= 3/2) = -0.0397 ^LeV, 

‘^E^AnJF = 5/2) = 0.0595 fieV. 

^^{j-i/2-^j-i/2),vp(^ — 1 / 2 ) — — 0.1111 lJ,eV, 

= 3/2) = -0.1757 ^2ev. 

They evidently show that two our approaches to the calculation of hyperhne structure in 
muonic deuterium P-states lead to the same results. Two-loop vacuum polarization correc¬ 
tions shown in Fig. [T] are calculated in a similar way. They are included in Appendix C. 
Their numerical value is essentially smaller (see Table |T|). 



a 




FIG. 2: Vacuum polarization effects in the second order perturbation theory. Dashed and wavy 
lines represent correspondingly the Coulomb and hyperfine interactions. 

For a completeness, we analyze vacuum polarization corrections of order in second 
order perturbation theory, which are determined by the reduced Coulomb Green’s function 
34 l35| (see the amplitude in Fig. [2](a)): 


G2p{r,r’) = 


IJ^^Za) / 3 
36z‘^z'‘^ IdTT 


nn 


' G-(2+^' 




(49) 


g{z, z') ^ 24 z^ + 36 sp> + 36 z^z^ + 24 zl + 362<2| + 362^2/. + 492 /. 2 > - Sz^z^- 
—12e^^(2 + 2< + 2p2^ — 32 p^ + 122^20[— 2 C + Ei(z^) — hiz^ — /n2>], 
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where C = 0.5772... is the Euler constant, z = hhr, 2 ;< = min{z, z'), z^ = maxjz, z'). Using 


flAlD and (H9|) . we obtain the following integral expression for VP correction [2l| 


dx 


oo g-o:' 


vp,soPT 24mim2 547r Jq Jq 


dx'e 


-x\ 1+ 


2meg 


X 


X 


TfT , ruiKd ^ ^ 

Ti + —7—— zTi - 1 + a^jTs 
7712(1 + Kd) 


(50) 


Similarly, the correction of vacuum polarization and quadrupole interaction in second order 
PT has the form: 


^^vp,Q,SOPT 


259277 



g{x,x')x 


(51) 


X 


(S 2 ■ S 2 ) - 3(s 2 ■ n)(s 2 ■ n) 


(dp i - \5p z + 5 ) • 0.1120 (peU), ] = f = \ 

(V^dp i - ■ 0.1120 (peU), ] = |, j' = i. 


The coordinate integration over x, x' is performed again analytically and numerically over 
Summary numerical values of contributions in the hrst and second orders of PT to the 
P-state energies are presented in Table [H Table [TTl separately for diagonal and off-diagonal 
matrix elements. 

Based on the amplitudes it is possible to find the nuclear structure correction 

(the index str designates this contribution) to hyperhne splittings. For this aim, we will 
introduce into them an additional factor (—r^k^/6) with the deuteron root mean square 
radius connected with the expansion of form factors and omit factors containing the deuteron 
magnetic moment. After evident simplihcations we obtain the following contributions to 
hyperhne splitting potentials for states 2 Pi /2 and 2 P 3 / 2 : 


zhfs 


Ti^,str{v, q)j=i/2(^ = 3/2; 1/2) = 


rhfs 


Zarl 

C"[P<1 

(pq)^ 

12 

7772 

pq 

Zarl 

/ \nQ- 

(pq) 

12 

^ 27772 

pq 


rhfs 


T'i7,sir(P,q)j=3/2(^ = 5/2; 3/2) = 


Zar^ r 57771 r (pq) 


]_2(£3)!+2(i+^)[p,+ (p3)!]}. 

pq 2 pq 

(52) 

pq 0 4 pq 

(53) 


-]-4 


"t t 


12 ^ 67772 ^ pq 


-pq]i 


20(pq)2 2, 


.(pq)' 


6 pq 3 4 pq 

(54) 


Further integration and consideration of the general normalization factor directly lead to 
the following splittings: 


= 1/2. F = 3/2; 1/2) = ZZll f !!!l + 


16777i7772 V^2 


= 0.0032 peV, 


AP/,^/,,,(j = 3/2,P = 3/2;l/2) = 


1 777i 


32r?7ir?72 V^2 


= 0.0016 /reU, 


(55) 

(56) 
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= 3/2. F = 5/2; 3/2) = - a„) = -0.0026 ,.V. (57) 

As expected, these corrections are very small and do not affect the comparison of theo¬ 
retical results and planned experimental data. Other corrections of order a® are discussed 
in two Appendixes B and C. 


TABLE II: Off-diagonal matrix elements in the hyperfine structure of P-wave muonic deuterium. 


Contribution to HFS 

2^Pi/2,3/2) 

2"^Pi/2,3/2) 

leading order 

correction 

-126.0372 

-199.2824 

quadrupole correction 
of order 

613.6320 

-194.0475 

vacuum polarization 
correction of order a® 

-0.1437 

-0.2271 

quadrupole and vacuum 
polarization correction 

of order 

0.0891 

-0.0282 

relativistic correction 

of order a® 

-0.0043 

-0.0067 

vacuum polarization 
correction of order a® 

0.0001 

0.0001 

Summary contribution 

487.5360 

-393.5918 


V. SUMMARY AND DISCUSSION 

In this work we investigate the hyperfine structure of energy levels related to the P-wave 
states of muonic deuterium on the basis of three dimensional quasipotential approach in 
quantum electrodynamics. To increase the accuracy of the calculation we take into account 
the leading order contribution and several basic corrections of order and a®. These cor¬ 
rections are connected with the vacuum polarization effect, quadrupole interaction, nuclear 
structure and relativistic effects. Some corrections are obtained in analytical form, but most 
part of contributions to the energy spectrum is presented first in integral form, and then cal¬ 
culated numerically. All results are presented in Tables IIIIITIIIIII giving the values of diagonal 
and off-diagonal matrix elements and the positions of the P-energy levels. 

We would like to point out three main our results obtained in this work. 

1. New approach based on the use of special type projection operators on the states 
with definite quantum numbers of atomic angular momentum F and total muon angular 
momentum j is developed. It allows to simplify essentially the construction of the particle 
interaction operator through the use of computer methods for calculating Feynman ampli¬ 
tudes. In particular, this method can be useful when working with different loop corrections. 

2. We have increased the accuracy of the calculation of P-wave hyperfine splittings 
primarily due to the correct account the corrections of the fifth order over a. To this end. 
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TABLE III: Hyperfine structure of P-states in muonic deuterium 


State 

Energy, meV , [161. 

Energy, meV 

2 ^Pi/2 

-1.4056 

-1.40534 

2"Pi/2 

0.6703 

0.67031 

22^3/2 

8.6194 

8.62002 

2"P3/2 

8.2560 

8.25618 

2®P3/2 

9.3729 

9.37183 


the contributions have been built into the operator of the interaction of particles that are 
connected to the vacuum polarization and quadrupole interactions. We check the obtained 
results in two ways, in the formulated framework of tensor projection operators in momentum 
representation and the more traditional method for the calculation of corrections in the 
energy spectrum in the coordinate representation. Moreover, in our calculation we take into 
account the contributions not only the hrst but also the second-order perturbation theory. 

3. New higher order 0{a^) corrections are calculated. These corrections although small 
numerically and do not affect on the comparison with future experimental data, but clarify 
the structure of the perturbation series for the hyperhne splittings. 

Let us present more detail comparison of the results with previous calculations in [l^ 22 


Being different in the method of obtaining corrections of leading order O(a^) our results 
coincide with 
0,122 


22| . We mean both the spin-orbit, spin-spin contributions of order 0{a 


and quadrupole corrections of the same order [16 


a corrections which are differ signihcantly from the results of [Ig. In paper [16 


But we obtain the hfth order in 
the vacuum 


polarization corrections to hyperhne part of the Breit Hamiltonian are determined by the 
following modihcation of the potential with Z > 0: 


IdV 
r dr 


Za 


1 + ^ y /3(?)<i?(l + 2m,ir)e 


( 68 ) 


This leads to appearing of special factor of the form (1 -|- e 2 p) with numerical value e 2 p = 
0.000391 for the quadrupole correction and 


^2P 


a 

37r 


" / 1 2az 

VO+TF (T+T> 


(59) 


for the Uehling correction to the Breit Hamiltonian. Numerically, the coefficient in (j5^ is 
equal to the same value e 2 p = 0.000391. In our calculation we demonstrate that the vacuum 
polarization corrections to P-states are determined by different potentials (compare (158|1 
with our formula (HSD) a nd have different form for states with various quantum numbers F 


and j. In contrast to [16(] we have performed exact construction of corresponding potentials 
for different P-states and obtained through them numerical results that can not be reduced 
to a factor fl5^ . Our results are checked by two independent methods. As an example, we 


give a comparison of our vacuum polarization plus quadrupole interaction contributions to 
hyperhne splittings of the level 2 P 3/2 with the results of |l^ . In [0 these contributions 


are equal to = 3/2; 1/2) = —3iJ,^a^Qd£2p/80 = —0.3058 /reH, AE^y^iF = 


hhfs 


0=3/2 


j=3/2\ 


5/2; 3/2) = fra‘^QdS2p/4:8 = 0.1699 fieV and diher essentially from our corresponding values 
(—0.4394) fieV and 0.2441 fieV. The same situation occurs for other VP corrections. 
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Summing all diagonal and off-diagonal matrix elements we obtain the following energy 
matrix 


2^P^ 


1/2 




1/2 


M = 


3/2 


2 ^P 


3/2 


2^P 


3/2 


2 -Pl/2 
-1381.5765 
0 

487.5360 

0 

0 


2‘^P 


1/2 

0 

690.7897 

0 

-393.5918 

0 


2 P 3/2 
487.5360 
0 

8596.2539 

0 

0 




312 

0 

-393.5918 

0 

8235.7070 

0 


2 ‘’ 7 ’ 3/2 

0 

0 

0 

0 

9371.8295 




(60) 


Its diagonalization leads directly to the position of the energy levels 2P (see Table m 
and hyperhne splitting intervals which can be measured in the experiment. Accounting the 
accuracy of the calculation, we have added one extra decimal place in our results in Table HTTl 
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Appendix A: Basic contributions to hyperfine structure in coordinate representation 


Basic contribution to hyperhne structure is determined by hyperhne part of the Breit 
Hamiltonian 


= ^“(1 + '‘^'1 [ 1 + + “4 , 

Zmim2r'^ m2\}-+ Kd) Zmim2r^ 

(Al) 

where mi, m 2 are the muon and deuteron masses, are the deuteron and muon anoma¬ 

lous magnetic moments, Si and S 2 are the spin operators of muon and deuteron, n = r/r. 
The operator (lAip does not commute with the muon total angular momentum J = L -|- Si. 
As a result there is the mixing between energy levels 2Pi/2 and 2 P 3 / 2 . 

For the calculation of diagonal matrix elements ^2Pi/2 2Pi/2^ and 


2P 


3/2 




2 P 3 / 2 / we use the Coulomb wave function of 2P-state in coordinate rep¬ 


resentation: 


^2p{r) = 0), W = ^Za. 


(A2) 


The angle averaging in flAll) can be carried out by means of the following replacements 


Si —V J -——, L —V J - 


J2 


J2 


(A3) 
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which give the eigenvalues of the corresponding operators: 




Hi +1) -;(; + !) + - 


(£'J)=2 


Hi + i) + i{i + i)-- 


i^SiJ STliTlj'j SLiLj SLjLi^. 


The diagonal matrix elements have the general form: 

t;+ 


J^hfs _ + Kd) 


48mim2 


miKd TfT , Yfr 

Ti — (1 + a^)T2 


^ 2(1 + Kd) 


(A4) 

(AS) 

(A6) 


where the operators Tj are defined in flT7)l . Substituting here Ti and T 2 for definite quantum 
numbers F and j, we obtain the leading order contributions to the hyperfine structure of 
2 Pi /2 and 2 P 3/2 states: 


2 771/1/5 _ 
-^1/2 “ 


aV^(l + Kd) 
\%m 1 m 2 

4 p/7/5 _ Q:V^(1 + Kd) 

36mim2 
aV^(l + i^d) 


-^1/2 “ 


2 p/7/5 _ 

-^3/2 - 


72mim2 
q:V^( 1 + /td) 
36mim2 

/ 7 /s _ + /trf) 

24mim2 


4 p/7/5 _ 
-^ 3/2 “ 


6 p77/5 _ 

-^3/2 “ 


1+; 

miKd 1 

2m2{l + Kd) 

~2 

= -1380.3360 fieV, 

(A7) 

1 + 

miKd 1 

2m2(l Kd) 

T_ 

= 690.1680 fieV, 

(AS) 

2 + 

5 miKd 
27712(1 -f Kd) 

Y 

= 8162.2889 fxeV, 

(A9) 

'2 

.5 + 

miKd 

27712(1 + Kd) 

10 

= 8583.2315 /ieV, 

(AlO) 

'2 

.5 + 

miKd 

27712(1 + Kd) 

10 

= 9284.8027 fxeV, 

(All) 


where we take into account the fine structure interval = 8.86386 meV calculated 


m 


16l. l37l |. All expressions flATp - flAlip contain the correction to the anomalous magnetic 
moment of the muon. 


Off-diagonal 
2Pi/2 I I 2P3/2> 


matrix 

F=3/2 


elements 


(2Pi/2 I I 2P3/2> 


F=l/2 


and 


are essential to achieve a high accuracy of the calcula¬ 


tion. They differ by the value of atomic angular momentum. The angular averaging by 
means of flASp leads to Ti = 2 T 2 . For the calculation (LS 2 ), we use the general formula for 
the matrix elements of the scalar product of two irreducible tensor operators: 

{fs 2 F I (pi ■ T2) I JS 2 F) = {-ir^’^'-^W{jS 2 fs 2 ; PI) {f II pi II j) (S 2 II p2 II S 2 > , (Al2) 


where W{js 2 j's 2 ] FI) is the Racah coefficient. Applying (1A12p to (LS 2 ) we find: 

{fs 2 F I (P ■ S 2 )\jS 2 F) = (_l)-.-F-5.+L+3/2+y + l)(2j + l)x 


X y/(2 s2 1)(s2 + 1)s2(2P 1)(P l)p| 1 } { j i 1 

Two off-diagonal matrix elements of the operator Pi have the form: 


1 1 
2’^’2 


[L ■ 82 ) 


3 1\ 

2’ ’2/ 3 ’ 


1 3 

2 ’ 2 


[L ■ 82 ] 


3 3\ 

2’ ’2/ 3 


(A13) 

(A14) 


V, (A15) 
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where the 6j-symbols are taken from 

Using flAlSp . we obtain the leading order contributions to off-diagonal matrix elements 
of the Breit Hamiltonian flAlj) : 


j^hfs,off—diag O /i (1 -|- ( \/2 

48mim2 1 6 


rphfs,off—diag CT /i (1 “1“ Kd'j / \/5\ 


1 + 


1 + 


2miKrf _ 
m2{l + Kd) 

2miKd 


= -126.0372 /ieU, (A16) 


= -199.2824 /ieU. (A17) 


m2{l + Kd) 

There exist higher order corrections to flA16p and flA17p which are related to additional 
interactions and examined above. 

Appendix B: Relativistic corrections to hyperfine structure 

Relativistic corrections of order a® can be calculated by means of the Dirac equation 


38l. l39[ . We present here only a sketch of the output of the hnal formula for the numerical 


estimate. In the Dirac theory the hyperhne part of relativistic Hamiltonian has the form: 

\r X cd 




egNlJ‘NS2- 


(Bl) 


where /iv is the nuclear magneton, gjsi is the deuteron gyromagnetic factor. To hnd the 
expectation value of fIBip over atomic wave functions we should use the Wigner-Eckart 
theorem expressing initial matrix element through the reduced matrix elements: 


r X a 


^W{js2fs2; FI) (sa || Sa || Sa) 

Calculating the hrst reduced matrix element we can simplify flB2p as follows: 


(B2) 


= egN^^Ni-l)^^^^'-W{‘2s2 + l)(sa + l)saV(2j' + l){f + l)jW(jsa/sa; FI) x 

\r X ckI 


xy 

In the case of diagonal matrix element we have: 
[r X a 




-1 




/■°° Ak 

jg, ) = -iAkkRkk, Rkk = 2 j gk(,r)fk{r)dr, - iA^k = ^^2 _i 


(B4) 

The radial matrix elements are calculated analytically with the use of exact Dirac radial 
wave functions. After their expansion over a we hnd |39|: 


mPin) = RCiPm) = 


+ (B6) 


24 


As a result, general expressions for relativistic corrections to diagonal matrix elements take 
the form: 

y X +1) - +1) - +1)1. (B6) 








































19 




a®(l + Kci)iJ? m\ 7 


48,n.™. (B7) 

Numerical results for separate P-states are presented in Table [H Relativistic corrections to 
off-diagonal matrix elements are evaluated in a similar way. The radial and angular integrals 
in this case take the form: 


^kk 


gdr)f3{r) + g3{r)fi{r) dr, - iAi = 


[(/ +1/2)2^ 


2 / + 1 


= V, (B8) 


where the indexes near radial wave functions designate the values of muon total angular 
momentum j. Radial integrations lead to analytical formulas and corresponding numerical 
results 


phfs,off-diag 

reZ,F=l/2 


phfs,off-diag 

reZ,F=3/2 


a®(l -|- ml 3\/2 
48mim2 32 

a®(l -|- 3\/5 

48mim2 32 


-0.0043 peR, 
-0.0067 fxeV. 


(B9) 

(BIO) 


Although their size is extremely small compared with other corrections we have included 
them in Table UTl bv inserting numerical values with an accuracy 0.0001 peR for the deh- 
niteness. It shows the relative numerical value of obtained corrections. 


Appendix C: Two-loop vacuum polarization corrections to hyperfine structure 


Two-loop vacuum polarization corrections presented in Fig. [T]b,c,d have the order a®. We 
divide them into two parts: loop after loop contribution (vp-vp) and two-loop contribution 
to polarization operator (2-loop vp). For their calculation we use corresponding potentials 
in coordinate representation constructed in the same way as in j2l|: 




1 + 


l'y,vp—vp V 

miKd 


O 3 \ ‘i I I / P(P)dPT2 -7 

2mim2r^ \onJ Ji Ji ^ — d 


X 


{L ■ S2)[e^(l + 2meer)e-2™=R - 7]\l + 2meVr)e-^^^^^]- (Cl) 


^ 2(1 + Kd). 

-(1 a^) ^4m2r2[^'^e"2™A7’ _ ^4g-2me»?rj ^ . n){s 2 ■ n)]-F 

-|-[,^^(1 -f- 2me^r)e“^”^'=‘^'’ — //^(l -|- 2mer]r)e~‘^^‘'^'^] x [(si ■ S 2 ) — 3(si ■ n)(s 2 ■ «-)] 


2 rl 


\yhfs / ^ ^ Za{l Kd) 2 Pa 

2-ioop vp 2 mim 2 r^ ^\'^J Jo ^ 


f(v)dv 

^ ^ e X 


X 


1 + 


miKd 


1 + 


2 mf,r 


[L-sd- 


-(1 -|- a^) 

+ U + 


7712(1 Kd) ) ' \/l - 1 

-^[(si ■ S 2 ) - (si ■ n)(s 2 ■ n)]-F 


1 — n 

27T7f.r 


\/l — 


[(si ■ S 2 ) - 3(si ■ n)(s2 ■ n)] 


(C2) 
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Averaging fICip and flC2p over the Coulomb wave functions we obtain their numerical values 
in the hyperhne structure which are presented in Table [I]-Table [TTl The muon vacuum polar¬ 
ization correction is evaluated by means of a replacement rUe —)■ rrii in (jl6]). Its numerical 
value also is included in Table HJ-Table UTl 

For the calculation of contributions in the second order PT we should use in basic ex¬ 
pression 

^EioPT,^ = 2 < . Q . ( 03 ) 


the potential corresponding to pure hyperhne interaction or to hyperhne inter¬ 

action corrected by the vacuum polarization ehect. Aa a second perturbation we use the 
Coulomb potential of one-loop or two-loop order. All resulting matrix elements are calcu¬ 
lated analytically in a standard way in the integration over the coordinates of the particles 
and numerically by spectral parameters. Other details of their calculation can be founded 
in our previous papers 0 0, iH. 

The two-loop vacuum polarization contribution to hyperhne structure of order a® is 
determined also by the third order PT. In this case we should use the following expression: 


AErpQpTvp — {'4^ n 


AV^P ■ G ■ AV^^^ ■ G ■ AV^p 


^n) + 


+ 2 

I AV^f^ I 'iPn) U 


AV^P -G-V^p-G- 


l4r. 


-2 (V’n I AV^p I Ipn) I'lpn 


AV^P -G-G- AV^p 
AV^P -G-G- AV’^f^ 


'ipr 


A 


(C4) 


Using further exact perturbation potential flAlD . modihcation of the Coulomb potential AU*" 
and the Coulomb Green’s function G fH9|) . we obtain numerical values of corresponding cor¬ 
rections which are written in Table [I] as a separate line. Numerically the vacuum polarization 
contributions of order a® are extremely small and will not have a signihcant impact on the 
comparison with future experimental data. 
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